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We investigate the possible molecules composed of two heavy flavor baryons such as "AqBq"(Q = b, c) 
within the one-pion-exchange model (OPE). Our results indicate that the long-range n exchange force is strong 
enough to form molecules such as [EgE' e ]£~y 2 (<2 = b, c), [LqAqY^J^Q = b, c), [£ 6 Ej]£lj and [H^S^]^ 
where the S-D mixing plays an important role. In contrast, the n exchange does not form the spin-singlet AqBq 
bound states. If we consider the heavier scalar and vector meson exchanges as well as the pion exchange, some 
loosely bound spin-singlet S-wave states appear while results of the spin-triplet AqBq system do not change 
significantly, which implies the pion exchange plays an dominant role in forming the spin-triplet molecules. 
Moreover, we perform an extensive coupled channel analysis of the AqAq system within the OPE and one- 
boson-exchange (OBE) framework and find that there exist loosely bound states of AqAq(Q = b, c) with 
quantum numbers I(J F ) = 0(0 + ), 0(CT) and 0(1~). The binding solutions of AqAq system mainly come from 
the coupled-channel effect in the flavor space. Besides the OPE force, the medium- and short-range attractive 
force also plays a significant role in the formation of the loosely bound A^A^ and AjA;, states. Once produced, 
they will be very stable because such a system decays via weak interaction with a very long lifetime around 
1(T 13 ~ l(T 12 s. 

PACS numbers: 12.39.Pn, 14.20.-c, 12.40.Yx 



I. INTRODUCTION 



In the past few years, many exotic charmonium-like states have been reported by the Belle, BARBAR, CDF and DO collab- 
orations, such as X(3872) fl, X(4160) 11,7(4260) H, and Z + (4430) H. Recently, the BELLE Collaboration observed two 
charged bottomonium-like resonances Zfe(10610) and Z;,(10650) in the hidden-bottom decay channels ^■ ± Y(n5') (« = 1, 2, 3) 
and i&hbimP) (m = 1, 2) of Y(5S) @|. 

Since many of these states do not fit into the conventional qq picture in the quark model 
easily, how to interpret these "exotic" states becomes a challenging problem. The prominent feature of these states is that they 
are near the threshold of two charmed or bottomed mesons. For example, X(3872) lies close to the threshold of D°D*° while 
Z/,(10610) and Z|,(10650) are near the threshold of BB* and B*B* respectively. Inspired by this striking feature of these exotic 
states, many physicists attempted to interpret them as hadronic molecules composed of heavy mesons. 

A hadronic molecular state is a loosely bound state of hadrons. Voloshin and Okun began to investigate the existence of the 
bound states composed of charmed meson and antimeson De Rujula et al proposed i/r(4040) might be a D*D* molecular 
state 10] . Tornqvist explored the possible deuteron-like meson-antimeson bound states with the pion exchange potential 
Liu et al. investigated the possible molecular states composed of heavy mesons within the framework of the one-boson-exchange 
model (OBE) HlOtl - Ding et al. also gave a dynamic study of meson-meson molecular states with the one-boson-exchange model 
at the quark level fllll . Sun et al. interpreted the newly observed Zfe(10610) and Z/,(10650) as BB* and B*B* molecular states 
respectively lfl2ll . 

Actually, the idea of the loosely bound molecular states is not new in nuclear physics. It's well-known that the deuteron is 
a very loosely bound state composed of a proton and neutron. The interaction between the proton and neutron comes from the 
color-singlet meson exchange. Besides the long-range attraction from the pion exchange, the S-D mixing, the medium-range 
attraction from the correlated two-pion exchange (or in the form of the sigma meson exchange), and the short-range interaction 
in terms of the vector meson exchange combine to form the loosely bound deuteron. 

It's quite natural to extend the same formalism to the heavy baryon sector. Since the heavy baryon contains a charm or bottom 
quark, its large mass reduces the kinetic energy and helps the formation of the bound states. Froemel et al investigated the bound 
states composed of heavy hyperon and nucleon by rescaling the nucleon-nucleon potential in lfl3ll . Julia-Diaz et al explored 
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the bound states composed of double-charmed hyperons ITill . In our previous work lfl9l . we performed a study of the systems 
A C A C (A C ), H C H C (H C ), £ C S C (£ C ), E' C E' C (Z' C ) and fi c Q f (Q c ). 

In the present work we shall study the systems with two different heavy flavor baryons. For simplicity, we denote the systems 
with two same baryons as "AgAg" and the systems with two different baryons as "Agfig". The difference between the two 
systems is that the "Agfig" system contains the contributions coming from the K*, K°, K°, K f± , K*° and K*° exchange while 
the "AgAg" system does not. 

Among the possible loosely bound states composed of a pair of heavy baryons, the AgAg (Q = b, c) system is particularly 
interesting since it is the heavy analogue of the well-known H dibaryon. Since it was proposed by Jaffe in Ref. |j5ll . there have 
been lots of theoretical and experimental efforts. Recent investigations include the Lattice QCD calculation IU6I1 . calculations 
using the chiral effective field theory lfl7ll and the quark model lfl8ll . In this work, we shall perform an extensive coupled channel 
analysis of the AgAg {Q = b, c) system and investigate the role of the OPE and sigma/omega/rho meson exchange in the 
formation of the possible loosely bound state. 

This work is organized as follows. After the introduction, we present the formalism in Section [II] which contains the La- 
grangians, the coupling constants and the effective interaction potentials. The formalism for the coupled channel analysis of 
the AgAg system is given in Section Hill In Sections ||V]and[Vj we show the numerical results for "AgBg" and AgAg systems, 
respectively. The last section |VI] is a brief summary. Some useful formulae and functions are given in the Appendix. As a 
byproduct, we also collect the numerical results for the loosely bound states composed of a pair of heavy baryon and anti-baryon 
in the Appendix. 



II. FORMALISM 
A. The Iagrangian 



The heavy flavor baryon contains a charm or bottom quark and a diquark (two light quarks). In the heavy quark limit 
(niQ — > oo), the charm or bottom quark can be viewed as a static color source. The SU(3) flavor symmetry of the baryon is 
determined by the diquark. The heavy flavor baryons can be classified in terms of the symmetry of the diquark. The symmetric 
one belongs to the 6-representation while the antisymmetric one belongs to the 3-representation. On the other hand, the spin of 
the diquark is either or 1 which is antisymmetric or symmetric under the exchange of its two light quark spins. The baryon 
is a fermion system. Its total wave function should be antisymmetric under the exchange of its two light quarks. Therefore, the 
spin and the flavor of the diquark are correlated with each other. Taking the color wave function into account, the diquark in the 
6-representation should be spin-triplet while the one in the 3-representation should be spin-singlet. The spin of the baryon in the 
6-representation is either | or | while the spin of the baryon in the 3-representation is only \. 

In the following, we follow the notations in Ref. |23i] and list the heavy flavor baryon matrices and the exchanged meson 
matrices. The heavy flavor baryons are 
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The spin | baryon is marked with *. The superscript is the third component of its isospin. Q = b or c denotes the corresponding 



heavy quark. The exchanged bosons are 
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The exchanged bosons include the pseudoscalar and vector mesons given in Eq. ([2]) and the scalar meson cr. The lagrangians 
built under the SU(3)-flavor symmetry read as 



where 



■C - XLphh + -Cvhh + -Ca-hh, 

-Cphh = 8 p b 6 bJi [B 6 iy 5 MB 6 j + g pBi B } Tr [B 5 iy 5 MB^ + {g pJ 5 6 B 3 Tr [B 6 iy 5 MB 5 ] + h.c] , 

£vhh = 8vb 6 bJt [^n] + ^Tr [Bea^VBe] + g.^Tr [s^Wj] + ^Tr [b^W* B 3 ] 



(3) 



3 



+ \g V B 6 B S Tr [B& M WB 5 ] + 2 ^fg_ Tr [5 6 <v5^2» 3 ] + A.c. j , (4) 

Tr [B 6 o-B 6 ] + go-^Tr [fljo-flj] . (5) 

A1 and 'V are the exchanged pseudoscalar and vector meson matrices, respectively, which are given in Eq. (0. B(, and B3 are 
the heavy baryon matrices shown in Eq. (HJ. me and tn% are the masses of the heavy baryons belonging to the 6-representation 
and 3-representation, respectively. "g P B 6 B 6 , 8vB 6 b 6 , ■ ■ ■" are the coupling constants. 



B. Coupling Constants 

We derive the coupling constants in Eqs. ([3]|5]l from those between the nucleon and the light meson within the Quark Model 
(QM). For the vector meson exchange, we adopt the lagrangian without the anomalous magnetic term at the quark level as in 
Ref. l24h . One can refer to Ref. lfl9ll for the specific expressions of the couplings at the quark level in terms of those between the 
nucleon and the exchanged mesons. We list the coupling constants we need below and collect their numerical results in TableU 

4 V2 nij +mt 2 V3 + m / 

pseudoscalar exchange, g pB( ,B 6 = 7 gnNN ~ , g P B 6 B s = =-g*ftw— = . 8 p b- 3 b s = 0. (6) 

5 zmjv ~ j Ivrifj 

2 1 

scalar exchange, ga-B 6 B 6 = ^ga-NN, ga-B,B } = -^go-NN- (7) 

r 4V2 V m ' OT / r 

vector exchange, g vB( ,B<, = 2S2g pNN , f vBBfs = ——{g pNN + f pNN ) 2y2g pNN , 

j niN 

gvB } B- } = ^2g p NN, fvB s B- } = - ^2g p NN, 

2 V3 V m ' m / 

gvB 6 B, = 0, JvB 6 B, = 7—(g P NN+f P NN) • (8) 

5 mjv 

In the above expressions, g„NN, go-NN, gpNN and f p ^ff are the coupling constants between the nucleon and the exchanged mesons. 
Their numerical values are known quite well. m# is the mass of the nucleon. m, and m / are the masses of the ingoing and outgoing 
baryons, respectively. From Eqs. (1618) , the value of the same coupling constant is slightly different for different systems if one 
takes into account the mass difference of the baryons of the same representation. For example, g P B 6 B 6 = ^-^gnNN "^"^ f° r me 

interaction vertex Y. c Y. c n while g P B 6 B 6 - ^-^gxNN — fm~ ' ^ or tne ^c^c 71 " vertex. The masses of the baryons and exchanged mesons 
are summarized in Table [TTI The values of the coupling constants are JUH^]: g% NN /4n = 13.6, g n NN = 0.4, g^. NN /4n - 5.69, 
£pAw/ 47r = °- 84 ' fpNNlg P NN = k p = 6.1, gi NN /4n = 20.0, fujNN I gwNN = K a = 0. In our case, we assume the SU(3) symmetry. 
Therefore we need only three couplings for the pseudoscalar, scalar and vector meson exchange respectively. We adopt the 
three independent ones as gnNN, ga-NN and g P NN, since they do not vary much among different models. Generally speaking, the 
physical results of the loosely bound deuteron system are very sensitive to the vector meson coupling constants. The recently 
proposed renormalization approach, which uses a regularized boundary condition, can decrease the dependence on the coupling 
constants Il27ll. 



C. Effective Potential 

Applying the lagrangians in Eqs.(|3]|5j one can derive the effective interaction potential in the momentum space. Given the 
hadrons are not fundamental particles, we employ a monopole form factor at each vertex to roughly describe the structure effect 
of the baryon 

A 2 - m? T A 2 - m 2 . 

f(O) = — = — (9) 

yW A 2 -Q 2 A 2 + Q 2 ' K> 

where A is the cutoff parameter by which we can regulate the exchanged momentum. m ex and Q are the mass and four momentum 
of the exchanged meson, respectively, and A 2 — A 2 - go- Making the Fourier transformation, 

W = f dQ 3 e iQr <V(Q)T 2 (Q), (10) 

{2ny J 
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TABLE I: The numerical results of the coupling constants, "x" means such a vertex does not exist. 
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one obtains the effective interaction potential in the coordinate space, which are given below. One can refer to Appendix IVIIAI 
for some formulae. Since the hadronic molecule is a loosely bound state, the hadrons are not expected to be very close to each 
other. We neglect the contact interaction piece 6(r) in the potential. The detailed information of the delta term is given in the 
Appendix. 

We expand the effective interaction potential in terms of — (m is the heavy baryon mass) up to order of . We also adopt the 
approximation 1 lm 2 A ~ l/m 2 B ~ l/(mAiriB) due to the large masses (ma, m B ) of the heavy baryons. The effective potentials have 
four terms: the central potential term "Vc, the spin-spin term 'Vss, the spin-orbit term "Vis and the tensor term 'Vj. 

Winn, a) = e V P s {r;n,a) + e V p (r;n,a), 
<V s (r,n=l,o-) = ( V s c (r;n=l,cr) + ( V s LS (r,n=l,cr), 

<V{nn,P) = r V' c (r,n,/3) + r Vl s (r,n,p)+ f V i ss (r,n,p) + r V > T (r, «,/?), 

where the superscripts, p, s and v, denote the pseudoscalar, scalar and vector meson exchange, respectively, a — n, rj, K*, K°, 
K° and f3 = to, p, <p, K**, K*°, K*°. n = 1, 2 denotes the direct and cross diagrams, respectively. For the scalar exchange, "n = 1" 
means that the cr exchange only occurs in the direct diagram. The specific expressions read 

pseudoscalar exchange, *y^(r;n,a;) = C p a — — — - — Ho(A, m a , r)cr A ■ (Tb, 

' 47r IzmAinB 

3 

<V p (r;n,a) = C p / l f 2p - " a H 3 (A, m a , r)S AB (r), (11) 
when u 2 = m 2 - (nif - mi) 2 < 0, they change into, 

<V p (r;n,a) = C p / l f 2p - °° M (A, m a , r)er A ■ cr B , 



<V p T (r;n,a) = C p J l f 2p ~° M 3 (A, m a , r)S AB (r), (12) 
47r UniAmB 
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TABLE II: The masses of the relevant mesons and baryons l28ll . The bottomed baryons S?, Hi and E' b (marked with *) have not been observed 
experimentally. Their masses are fixed to be: m E o = [m^+ + m s -]/2 and m-> = in-'- = [m^ b + m n J/2. 
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FIG. 1 : The exchanged mesons for different systems. 
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with 6 2 a = - \ml - (nif - ra,) 2 ], 
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C p n a , O nB and C s n ^ are the isospin factors. Their numerical values are given in Table HTT1 and the exchanged mesons are shown in 
Fig.Q] L is the relative orbit momentum operator between the two baryons "Ag" and "fig". Sa, Sb are the spin operators of the 
two baryons while S = Sa + Sb is the total spin operator. S ab{?) = 3oa ■ rer B • r / r 2 - <Ta ■ <Tb is the tensor operator. g\ p , g 2p , . , . 
are the coupling constants given in Eqs. dSJSJ. The values of m, read 

direct diagram, u a = u a = u B = 0, 

9 9 9 9 T 9 9l 9 9 9 

cross diagram, u a = m a - {m A - m B ) , 8" a = - \m a - (m A - m B )~\ , u B =m B - {m A - m B ) . (15) 



Substituting the masses of the corresponding baryons for m A and mj in Eqs. ( fTTTTBl i. one obtains the effective interaction 
potentials. Besides Eqs. (II 1 HI 2b in the one-pion-exchange (OPE) model, we need include the contributions from the other 
heavier exchanged mesons in the the one-boson-exchange (OBE) model. The potential within the OBE model reads 
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The systems with two spin-half particles are either spin-singlet (S=0) or spin-triplet (S=l). For the spin-singlet, we focus on 
the ground state 'So while for the spin-triplet, we take both 3 S i and 3 D; into account. The wave functions can be expressed as 



= y%(r)\ 1 S >, ¥(r) r = ( y « (r) )| 3 Si>+( 







j Dt >, 



where y s s (r), y$(r) and y T D are the radial wave functions. The operators can be written in the following matrix form: 
spin-singlet, cr A ■ o" B = -3, L ■ S = 0, L ■ Sa = 0, L ■ Sb = 0, S A b(?) = 0, 



(17) 



/io\_.„/oV8\,_/oo\, c (o \ , c (0 \ 

spin-triplet, cr A • cr B = I Q { J, S AB (r) = I ^ _^ I » L ■ S = I Q _ 3 J, L ■ S A = y Q _3 J, B ~ _2 J ■ ( 18 ) 



III. THE COUPLED CHANNEL ANALYSIS OF THE A G A e SYSTEM 

The AqAq [0(0 + )] system is very interesting, which can be viewed as the heavy analogue of the H dibaryon. The heavy 
quark mass m& jC , the S-D wave mixing and the coupled channel effect in the flavor space all may play an important role in the 
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TABLE III: The isospin factors. The superscript in the first column is the isospin of the system. The values outside the square brackets are for 
the direct diagram while the ones inside the square brackets are for the crossed diagram. 
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formation of the possible loosely bound states. Investigation and comparison of the A^A/,, A e A c , and AA systems may reveal 
which underlying mechanism is dominant. 

In the present work we shall perform an extensive analysis of AqAq with quantum numbers I(J P ) = 0(0 + ), 0(0~) and 0(1). 
We list the flavor channels which we take into account in Table II VI Besides the lagrangians given in Eqs. we also need the 

following effective Lagrangians: 



pseudoscalar exchange, X. P = 

scalar exchange, £ s 
vector exchange, 



Si 
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' Tr(5^ ryT5 5 v MB^) 
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+ {-i V2Ajrg v Tr [b^ (3" V v - d v V) y v y 5 B 3 ] + h.c.} 



V6 



(19) 
(20) 
(21) 



(22) 



The coupling constants are g 4 = 0.999, g 3 = V6g 4 , g 5 = - V2g 4 , ^ = 6.2, = 92.3 MeV, (B s g v ) = 12, (/l s g v ) = 19.2 GeV 1 
and (Aig v ) = —(A s g v )/ V8 If20ll . Besides the potentials in Eqs. (II 1II14I) . we also need the following potentials 



<V p {r) = C"(i, j) — [H 3 (A, m ex , r)A ten + H (A, m ex , r)A ss ] , 
247r 

3 

V(r) = C\(i, j) —H (A, m ex , r) + C 2 (i, j) — [-H 3 (A, m ex , r)A ten + 2H (A, m ex , r)A ss ] , 

47T YL7Z 
Til 

<V '(r) = C s (iJ)-^Ho(A,m IT , r ). 
An 



(23) 

(24) 
(25) 



where A len and Ass denote the tensor and spin-spin operators respectively. They are channel-dependent. Their specific expres- 
sions are given in Table M 

For the baryon masses, we use = 2518.0 MeV and mj? = 5832.5 MeV Ir28tl . Due to the conservation of the energy and 
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momentum, we keep the non-vanishing zeroth component of the exchanged four momentum Qq and define u as the following 



S e 2g <-> TqLq, u 2 = m 2 ex 






, AqAq <-> I g S* g , m 2 = m 2 A . - 




2 


- 










mi, - mi 


2 


/ 2 2 > 

mi:, - mi 1 


2 


EqIq <-» SqSq, m 2 = m 2 ^ 




, S* G Z* e <-> I e Z* G , m 2 = m 2 x - 


, Other channels, u 2 = m 2 x . 






V h Q J 


TABLE IV: The flavor channels for the AgAg system, where Q = b or c. 


Channels 1 




2 


3 4 


5 


6 7 


= 0(0 + ) A e A e ( 1 5 ) 




SgSg('5o) EgSg( 5 Z? ) 


J£r e ( s Z) ) 




/(J p ) = 0(0") AqAq( 3 P ) 










I(J P ) = 0(1") AflAcfPi) 


S e 2 e ( 3 Pi) 




2 e 2* e ( 3 / > i) 





IV. NUMERICAL RESULTS FOR THE A Q B Q SYSTEMS 

In our numerical analysis, we apply the Fortran program FESSDE l29ll to solve the multichannel Schrodinger equation. 
Solving the Schodinger equations with the potentials derived in the previous sections, we obtain the numerical results including 
the binding energy (B.E.), the root-mean-square radius (r rms ) and the probabilities of the individual channels. We also plot the 
dependence of the binding energy on the cutoff parameter in the Appendix. 

The hadronic molecule is a loosely bound state. Its constituents are expected to be well separated. One expects that the size 
of the molecules should be much larger than that of the conventional qq and qqq hadrons. Recall that the size of the deuteron 
is about 1.96 fm IT22tl . Tornqvist argued that the size of the meson-meson molecule is even up to 3 fm J2lll . We expect the 
size of the hadronic molecules composed of two heavy baryons should be comparable to the size of the deuteron. The size of 
the molecular system may tell us whether the present framework and numerical results are self-consistent or not. To be more 
specific, the size of the molecular states composed of two charmed (or bottomed) baryons is expected to be larger than that of 
J /if/ (or Upsilon). 

Generally the value of the cutoff parameter is determined through fit to experimental data. In our case, there is almost no 
information on the heavy baryon-baryon interaction through which we can extract the cutoff parameter. Fortunately, the one- 
boson-exchange model is very successful in explaining the deuteron with the cutoff parameter 0.80 GeV< cutoff <1.5 GeV, 
which provides us a good benchmark. 

The present OBE model is rather crude. For example, we adopt the same cutoff parameter for all the meson exchange. We 
plot the interaction potentials of the deuteron with the OBE model in Fig. [2] We present the numerical results in Table IVTl 



TABLE V: The specific expressions of operators A 1e „ and Ass f° r the individual channels. 5, and |cr rs = -jS tll aS m are the transition matrix 
and the spin operator of the spin-| baryons, respectively. One can refer to Ref. i20ll for their definitions. 



Channels A K „ Ass 



AqAq <-» XgSg 


3S* 


■K 


f - 


S f 






S l2 


A Q A Q <-> ZqL* q 


3<T[ 




f - 


(T\ 




f 1 


S <2 








f - 




• 0"«2 




• CT rs2 




3<xi • 


rcr rs2 


f - 






0"i • 


&TS2 




3S* 


•fS, 2 


f - 




S,2 




■S e 
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2 4 6 8 10 12 14 2 4 6 8 10 12 14 2 4 6 8 10 12 14 

rtGeV" 1 ] r[GcV~'] rfGeV" 1 ] 

FIG. 2: The interaction potentials of the deuteron within the OBE model. V\\, Vn and V22 are for the transitions 3 5i<-» 3 5i, 3 Si<-> 3 Di and 
3 D 1 <-» 3 D,, respectively. 



TABLE VI: The numerical result of the deuteron with the OBE potential. "A" is the cutoff parameter. "B.E." means the binding energy. "Pg" 
and "Pd" are the probabilities of the S wave and D wave respectively. 





without contact term 








with contact term 






A (GeV) 


B.E. (MeV) 


r rms (fm) 


Ps(%) 


Pd(%) 


A (GeV) 


B.E. (MeV) 


rrms (fm) 


P S (%) 


Pd(%) 


0.80 


9.31 


2.14 


93.19 


6.81 


0.80 


1.87 


4.06 


95.00 


5.00 


0.85 


18.77 


1.65 


92.40 


7.60 


0.85 


2.58 


3.59 


94.14 


5.86 


0.90 


29.45 


1.39 


92.14 


7.86 


0.90 


2.88 


3.37 


93.82 


6.18 



In order to study the effect of the contact interaction, we compare the results (1) when the 6(r) function is omitted and (2) 
when the S(r) function is explicitly kept. Without the contact interaction piece, the binding energy of the deuteron is 9.31 
MeV and the root-mean-square radius is 2.14 fm when the cutoff is 0.8 GeV. If we include the contact interaction, the binding 
energy decreases to 1.87 MeV and the root-mean-square radius increases to 4.06 fm. In other words, both approaches roughly 
reproduce the qualitative feature of the loosely bound deuteron. In the following, we present the numerical results without the 
contact interaction. For comparison, we collect the results with the contact interaction in Appendix I VII B I 

However, if we shut down the 3 D\ channel, we can not find the binding solutions, which means that the S-D mixing effect is 
very important in the formation of the loosely deuteron bound state although the probability of the D wave is as small as ~ 6%. 
One can refer to Fig. |3]for the variations of the binding energy and the root-mean-square radius of the deuteron with the cutoff 
parameter. 

It's interesting to investigate whether the long-range pion-exchange interaction plays a dominant role in forming the hadronic 
molecules. Therefore, in the first part we give the numerical results with the one-pion-exchange (OPE) potential for the systems 
where the pion exchange is allowed. In the second part we take into account the scalar and vector boson exchanges, which 
account for the medium- and short-range interactions, as well as the n exchange. 




A ( GeV) A (GeV) 

FIG. 3: The variations of the binding energy and the root-mean-square radius of the deuteron with the cutoff parameter. 
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TABLE VII: The numerical results for the spin-triplet system (S=l) with the OPE potential, "x" indicates no binding solutions. A is the cutoff 
parameter. B.E. is the binding energy while r nns is the root-mean-square radius. P$ and P D indicate the the probabilities of the S wave and the 
D-wave, respectively. Q = c or b denotes the charmed or the bottomed systems. 



S=l 

Q=c Q=b 



Systems 


A(GeV) 


B.E.(MeV) 


r rms (fm) 




Pd{%) 


A(GeV) 


B.E.(MeV) 


r rms (fm) 


Ps(%) 


Pd(%) 


ry ,7=1/2 


1.20 
1.30 
1.50 


0.09 
0.84 
5.15 


5.99 
3.51 
1.67 


97.82 
95.69 
92.15 


2.18 
4.31 
7.85 


0.90 
1.10 
1.30 


4.28 
13.63 
31.33 


1.53 
1.01 
0.74 


88.01 
85.20 
83.37 


11.99 
14.80 
16.63 


[2 e E e ]^," 3/ - 


2.00 
2.20 
2.40 


0.21 
4.31 
14.85 


5.07 
1.61 
0.98 


96.66 
89.98 
85.41 


3.34 
10.02 
14.59 


1.10 
1.20 
1.30 


0.96 
2.99 
6.44 


2.34 
1.53 
1.17 


86.98 
82.01 
78.41 


13.02 
17.99 
21.59 




1.30 
1.40 
1.50 


0.55 
1.76 
3.99 


4.20 
2.71 
1.94 


97.46 
89.71 
87.29 


7.54 
10.29 
12.71 


0.90 
1.00 
1.10 


8.61 
13.95 
21.54 


1.40 
1.16 
0.98 


71.76 
70.04 
68.35 


28.24 
29.96 
31.65 




1.80 
1.90 
2.00 


0.88 
2.59 
6.35 


3.25 
1.96 
1.42 


93.65 
90.33 
87.63 


6.35 
9.67 
12.38 


0.90 
1.10 
1.30 


1.38 
5.98 
16.55 


2.37 
1.35 
0.93 


83.17 
77.00 
72.40 


16.83 
23.00 
27.60 


[E Q E Q ]'f l 






X 






2.20 
2.30 
2.40 


0.46 
1.20 
2.39 


3.01 
1.99 
1.48 


94.06 
91.76 
89.69 


5.94 
8.24 
10.31 



A. The Results of The "A Q B Q " Systems with The OPE Potential 

From Fig. Q] and Table [TTIJ we can see that there exists the one-pion-exchange force for five systems: ['LqE'q] < -^ 1 ^ 2 ' 3 ^ 2 \ 

[2gAg]^ =1 and ["EqE!q]^°' 1 \ For [SgSg]^ 1 ^ 2 and [EgHg] 3 ^ 2 , the pion exchange exists in the direct channel while the pion 
exchange occurs in the cross channel for the other three systems. 

There are no binding solutions for the five spin-singlet (S=0) systems. For the spin-triplet (S=l) case, we list the numerical 
results in Table lVUl and plot the dependence of the binding energy on the cutoff parameter in Fig. [5] We obtain binding solutions 
for all the states except for [HcH^,] 7 ^, see Table IVIll In the charmed sector, loosely bound states [1, C 'E' C ] I S ^( 2 and [SfAc]^ have 
small binding energy around a few MeV for a reasonable cutoff parameter about 1.20 ~ 1.50 GeV. They are good candidates of 
molecules. For these two states, the D wave contribution is less than 13%. For the states [L C 'E' C ] , S '^ 2 and [EcE^,] 7 ^, the binding 
solutions exist with the cutoff parameter around 2.0 GeV and 1.80 GeV, respectively. And, the D wave probabilities for the two 
states are less than 15%. 

The bottomed case is similar to the charmed case except that the binding energy of the bottomed bound states is deeper. This is 
mainly because that the larger mass of the bottomed baryon reduces the kinetic energy. For the states [S^S^ 7 ~J /2 ' 3 ^ 2) , [ZbAt,] 1 ^ 
and [EfeSJ 7 ^ , we obtain loosely bound states with binding energy less than 30 MeV and the root-mean-square radius larger than 
0.7 fm, when the cutoff parameter is about 0.9 ~ 1 .30 GeV. The D wave contribution is larger for the bottomed systems than for 
the charmed systems, see Table lVIil In [30], the authors performed a study of the 1,qAq system at the quark level and obtained 
larger binding energy. The [E^E^jtl'j bound state appears when the cutoff parameter is 2.20 GeV. When 2.20 GeV< cutoff <2.40 
GeV, the binding energy is 0.46 ~ 2.39 MeV. 

By comparing the numerical results of the two pairs of isospin multiplets, [EgHg] 7-1 ' 2 and [SqH'^] 7-3 ' 2 and [SgH^]^ and 
[EqEq]^ 1 , we can see that the results are different for different isospin multiplets of the same flavor system since the potentials 
are isospin-dependent. Comparing the results of the charmed systems with those of the bottomed systems, it is obvious that the 
large heavy quark mass is salient in the formation of the molecular states. 
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B. The Results of The "A e B e " Systems with The OBE Potential 

In the previous subsection, we give the numerical results with the n exchange potential which accounts for the long-range 
interaction. Actually, only five systems out of the thirteen ones allow the n exchange. We find that the n exchange is not strong 
enough to form bound states for all the five spin-singlet systems. In order to make the individual role of the exchanged boson 
clear, we also give the numerical results within the one-boson-exchange model (OBE) in Tables lVIIlllIXl and plot the dependence 
of binding energy on the cutoff parameter in Figs. (|6]|9]l. 

' 7=1/2 

In the spin-singlet case, we still find no binding solutions for the state [EgH e ] s=( j (Q=b,c) even if we add the contributions 

' 7=1/2 

of the heavier vector and scalar meson exchange. Therefore, our results disfavor the existence of the molecules [ZgHg] J=0 
(Q=b,c). However, for the other isospin multiplet of this state with / = |, we find binding solutions for both the charmed and 
bottomed cases. A bound state of [S C H C ] S=1 appears with binding energy about 3.54 ~ 67.46 MeV when the cutoff parameter 
is around 0.8 ~ 1.0 GeV. The binding energy of the corresponding bottomed state [E^E' ]£l 3 ' 2 is about 156.78 MeV when the 
cutoff parameter is 1 .0 GeV. Such a large binding energy seems too deep for a loosely bound molecular state. For the [SgAg]^, 
system, the binding energy of the state [EcAJ^L 1 ,, is 0.28 ~ 47.34 MeV with the cutoff parameter 0.90 -1.10 MeV while the 
binding energy of ["EbAbYf^ is 0.34 ~ 62. 1 3 MeV with a cutoff parameter 0.80 ~ 1 .00 GeV. We obtain bound states for both the 
charmed and the bottomed cases for [S^Eg] with 7 = and 7=1. 

7=1 /2 

For the other systems without the n exchange, we also find binding solutions. The most interesting one may be [EqAq] s _ , 
which allows the cr and to exchanges in the direct channel and K* /K* exchange in the cross channel. For the charmed case, a 
very loosely bound states with binding energy 1.91 ~ 3.03 MeV appears when the cutoff parameter is 1.10 ~ 1.50 GeV. For the 
bottomed case, a bound state emerges with binding energy 10.33 ~ 28.65 MeV when the cutoff parameter is between 0.90 GeV 
and 1.50 GeV. They are very good molecule candidates. For the states [LcScYf^ 2 and [Xcfic]^ 1 ,,, we also obtain small binding 
energies and large root-mean-square radii with reasonable cutoff parameter 1 .0 ~ 1 .50 GeV as shown in Table IVIIII Our results 

' I— 1 /2 

are in favor of the existences of these molecular states. The binding energy of [S b Cl b Y s ~^ is 80.49 ~ 107.01 MeV with cutoff 
parameter 0.90 ~ 1.00 GeV. Again, such a large binding energy seems too deep for a loosely bound molecular state. 

In the spin-triplet sector, it is interesting to compare with the deuteron case. We plot the interaction potential of the 
£ e S^, ^(.Z^) = 5(1 + )] system in Fig. [4] From Fig. [2] and Fig. [4] it is clear that the potentials of the two systems are similar. 

Their binding solutions are also similar except that the E C E C system has even shallower binding energy and smaller D wave 
probability, as can be seen from Tables IVT1 and IIXI 

' 7=3/2 

The bound state of [E C E C ] S=1 disappears if we take the heavier scalar and vector meson exchanges into account. There is still 
no binding solution for the state [S C B C ]^ = L 1 1 when we consider all the contributions of the exchanged mesons. The binding energy 
of the state [LqAq]'^ becomes shallower in the OBE model. For the other systems with the n exchange, the numerical results 
within the OBE model are similar to those within the OPE model except that the binding energy becomes deeper as shown in 
TableHS 

From Table HXl one can see that there is no S-D mixing for the two states [EgAglj" 1 / 2 and [AqQ.q] 1 ^. Actually, for these 
two systems the results are the same for both the spin-singlet and spin-triplet cases because the potential is the same. For the 

' 7=1/2 7=3/2 ' 7=1/2 

states [H c i2 c J„_ 1 ' and [E C E C ] S=1 , there are no binding solutions. However, a very loosely bound state [E c A c ] s= j exists with 
binding energy 0.17 ~ 0.69 MeV when the cutoff parameter is around 1.00 ~ 1.40 GeV. The binding energy of its bottomed 

' 7=1/2 

counterpart [a b A},\ s _{ is 17.64 ~ 23.91 MeV with cutoff parameter around 1.0 ~ 1.40 GeV. We also obtain a loosely bound 

7=1/2 

state [H c n c ]„_j with binding energy 2.63 ~ 4.51 MeV when the cutoff parameter is between 1.00 GeV and 1.20 GeV. Once 
these three molecule states are produced, they should be very stable because their constituents Aq, Eg, Eg and Qg decay via 
weak interaction. Comparing the results of the OBE model with those of the OPE model, one can see that the contribution of 
the D wave decreases if we take into account the scalar and vector meson exchange, which implies that the S-D mixing mainly 
comes from the n exchange. 

V. NUMERICAL RESULTS FOR THE A Q A Q SYSTEM 

We investigated the AqA q \l {J P ) = 0(0 + )J system with the cr and a> exchange potential, but without the coupled-channel effect 
in the flavor space in Ref. lU9"tl . We find no binding solutions. Later, the authors in Ref. rf20ll considered the coupled-channel 
effect and studied this state using the pion-exchange potential. They found a bound state solution. It is intriguing to study the 
variation of the bound state solution with the heavy quark mass, the S-D mixing effect, the long-range OPE force and medium- 
/short-range interaction respectively. In the first subsection, we shall present the numerical results for the OPE model with the 
coupled-channel effect. In the second subsection, we will add the scalar and vector meson exchange force which also contributes 
to the transition in the flavor space. 
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TABLE VIII: The numerical results for the spin-singlet (S=0) case with the OBE potential, x means no binding solutions exist. 

S=0 

Q=c Q=b 
Systems A(GeV) B.E.(MeV) r„,„(fm) A(GeV) B.E.(MeV) r„,, t (fm) 



rv 77' i'='/ 2 



ry -y' ,7=3/2 



0.80 
0.90 
1.00 



3.54 
14.53 
67.46 



1.97 
1.14 
0.67 



0.80 
0.90 
1.00 



25.63 
53.38 
156.78 



0.74 
0.57 
0.40 



0.90 
0.95 
1.00 



28.50 
29.65 
44.09 



0.88 
0.86 
0.75 



0.90 
0.95 
1.00 



80.49 
80.91 
107.01 



0.49 
0.49 
0.45 



[~Lq£Iq\' s 1 



1.00 
1.20 
1.40 



0.56 
13.26 
36.77 



3.93 
1.09 
0.75 



0.90 
1.00 
1.10 



5.12 
18.07 
36.67 



1.19 
0.77 
0.60 



[s e A e r /2 



1.10 
1.30 
1.50 



1.91 
3.03 
2.87 



2.47 
2.03 
2.08 



0.90 
1.20 
1.50 



10.33 
28.65 
27.78 



0.96 
0.67 
0.68 



0.90 
1.00 
1.10 



0.28 
14.81 
47.34 



5.48 
1.09 
0.73 



0.80 
0.90 
1.00 



0.34 
16.02 
62.13 



5.19 
0.79 
0.52 



1.00 
1.30 
1.50 



4.13 
20.99 
26.92 



1.77 
0.95 
0.86 



0.90 
1.00 
1.10 



11.85 
32.21 
49.64 



0.90 
0.64 
0.55 



0.90 
1.00 
1.10 



1.29 
9.15 
33.32 



2.91 
1.32 
0.82 



0.90 
0.95 
1.00 



19.88 
29.69 
46.46 



0.78 
0.68 
0.58 



[s e A C r /2 



0.90 
1.00 
1.10 



0.58 
7.08 
24.23 



3.99 
1.47 
0.93 



0.90 
0.95 
1.00 



16.62 
26.45 
40.88 



0.83 
0.71 
0.61 



0.95 
1.00 
1.05 



6.67 
23.80 
44.48 



1.48 
0.92 
0.74 



0.90 
0.94 
1.00 



5.67 
28.80 
73.68 



1.19 
0.67 
0.49 



[3 3 e ]£ 1 



0.90 
1.00 
1.10 



8.18 
22.78 
56.04 



1.38 
0.95 
0.69 



0.90 
0.95 
1.00 



40.88 
53.05 
73.19 



0.61 
0.56 
0.50 



[E q Qq]' s 1/2 



0.90 
0.95 
1.00 



2.19 
12.13 
30.45 



2.31 
1.19 
0.86 



0.80 
0.90 
0.94 



1.44 
20.88 
40.58 



2.06 
0.78 
0.62 



[A e n e ]f° 



1.00 
1.10 
1.20 



5.40 
16.55 
32.02 



1.60 
1.04 
0.82 



0.90 
1.00 
1.10 



14.10 
36.07 
64.89 



0.86 
0.62 
0.51 
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TABLE IX: The numerical results for the spin-triplet (S = l) case with the OBE potential, x indicates no binding solutions exist. 

S=l 



s 


A(GeV) 


B.E.(MeV) 


Q=c 
r rmj (fm) 


Ps(%) 


Pd{%) 


A(GeV) 


B.E.(MeV) 


Q=b 

r„„.,(fm) 


Ps(%) 


Pd{%) 




0.80 


2.59 


2.28 


96.48 


3.52 


0.80 


22.96 


0.85 


90.73 


9.27 


1/2 


0.90 


15.15 


1.17 


95.02 


4.98 


0.85 


34.88 


0.74 


90.16 


9.84 




1.00 


55.44 


0.74 


95.17 


4.83 


0.90 


56.02 


0.63 


90.37 


9.63 














0.80 


0.46 


3.22 


90.74 


9.26 


3/2 






X 






0.90 


10.34 


1.08 


82.64 


17.36 














1.00 


1.04 


2.40 


86.31 


13.69 














0.90 


0.88 


2.38 


96.64 


3.36 


=1/2 






X 






1.00 


23.30 


0.76 


89.72 


10.28 














1.10 


10.78 


0.98 


89.47 


10.53 




0.90 


1.20 


2.93 


99.91 


0.09 


0.80 


3.07 


1.47 


99.73 


0.27 


= 1 


1.00 


10.04 


1.25 


99.67 


0.33 


0.90 


19.30 


0.78 


99.24 


0.76 




1.10 


25.78 


0.89 


99.39 


0.61 


1.00 


46.47 


0.58 


98.71 


1.29 




1.00 


0.72 


3.69 


100.00 


0.00 


0.90 


10.33 


0.96 


100.00 


0.00 


=1/2 


1.30 


3.03 


2.03 


100.00 


0.00 


1.30 


29.55 


0.66 


100.00 


0.00 




1.50 


2.87 


2.08 


100.00 


0.00 


1.50 


27.78 


0.68 


100.00 


0.00 




0.80 


0.01 


6.58 


96.76 


3.24 


0.90 


19.85 


0.98 


80.86 


19.14 


= 1 


0.90 


0.46 


4.49 


95.26 


4.74 


1.10 


9.56 


1.36 


77.43 


22.57 




1.00 


0.03 


6.48 


96.59 


3.41 


1.30 


3.35 


2.12 


73.50 


26.50 




0.90 


1.63 


2.62 


99.95 


0.05 


0.80 


3.78 


1.37 


99.90 


0.10 


1/2 


1.00 


17.30 


1.04 


99.88 


0.12 


0.90 


21.96 


0.75 


99.60 


0.40 




1.10 


49.98 


0.72 


99.90 


0.10 


1.00 


62.92 


0.53 


99.59 


0.41 














0.90 


5.45 


1.19 


97.78 


2.22 


3/2 






X 






1.00 


5.52 


1.19 


96.60 


3.10 














1.10 


1.31 


2.06 


97.61 


2.09 




1.00 


0.17 


5.44 


99.94 


0.06 


1.00 


17.64 


0.80 


98.81 


1.19 


=1/2 


1.20 


0.69 


3.71 


99.93 


0.07 


1.20 


23.91 


0.71 


99.41 


0.59 




1.40 


0.36 


4.60 


99.97 


0.03 


1.40 


23.14 


0.71 


99.92 


0.08 




0.90 


3.10 


2.08 


97.44 


2.56 


0.80 


11.16 


1.06 


88.45 


11.55 


'==0 
J 


1.00 


13.75 


1.16 


98.60 


1.40 


0.90 


24.87 


0.77 


94.63 


5.37 




1.10 


32.31 


0.84 


99.38 


0.62 


1.00 


51.15 


0.58 


98.27 


1.73 














0.90 


3.40 


1.47 


96.48 


3.52 


=1 






X 






1.00 


13.70 


0.89 


95.83 


4.17 














1.20 


8.11 


1.08 


96.54 


3.46 




1.00 


2.63 


2.09 


99.37 


0.63 


1.00 


28.78 


0.69 


95.78 


4.22 


= 1/2 


1.10 


4.51 


1.08 


99.33 


0.67 


1.10 


35.44 


0.63 


96.52 


3.48 




1.20 


2.73 


2.05 


99.57 


0.43 


1.20 


29.48 


0.67 


97.72 


2.28 




1.00 


5.40 


1.60 


100.00 


0.00 


0.80 


1.62 


1.88 


100.00 


0.00 


=0 


1.10 


16.54 


1.04 


100.00 


0.00 


0.90 


14.10 


0.86 


100.00 


0.00 




1.20 


32.02 


0.82 


100.00 


0.00 


1.00 


36.07 


0.62 


100.00 


0.00 
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. A=0.S0GcV 
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- A=0.90GeV 



10 12 14 



r[CeV-'] 



FIG. 4: The potential of system ~L C Z' C [l(J F ) = within the OB E model. Vn, Vn and V 22 are for the transitions 3 Si*^> 3 S u 3 Si^ 3 £>i and 

^Di^Df, respectively. 

A. Numerical Results for The A e A e System with The OPE Potential 

The state AgAg(Q=b,c) is the heavy analogue of the H dibaryon. We use the n exchange potential and include the coupled- 
channel effect. Actually, the pion exchange is forbidden for the AgAg system due to the isospin conservation of the strong 
interaction. The binding solution is mainly due to the coupled-channel effect. We give the numerical results in Table IXl and plot 
the dependence of the binding energy on the cutoff parameter in Fig. [TU] 

AgAg [l(J p ) = 0(0 + )]. We reproduce the numerical results of Ref. |j20| for the state A c A c ( l S Q ), and list them in Table [X] 

We also extend the same formalism to the bottomed sector. A bound state A^A^'So) with binding energy 9.09 ~ 26.99 MeV 
appears when the cutoff parameter is chosen between 0.80 GeV and 0.90 GeV. Correspondingly, the root-mean-square radius 
varies from 1.08 fm to 0.77 fm. However, if the cutoff parameter increases to 1.10 GeV, the binding energy will reach as high as 
105. 17 MeV. There are five channels shown in Table HVIfor this state. The A/jA^'So) component is dominant with a probability 
about 90%. The contribution of the components 'Eb'^iA 1 S o) and Z^Zj^'So) is very small, around 1%. 

AgAg |/(./ p ) = 0(0 )J . We also extend the same analysis to the case with the orbital excitation L = 1 and obtain a loosely 
bound state of A f A c [0(0~)] with binding energy 0.19 ~ 24.97 MeV when the cutoff parameter is chosen between 1.36 GeV 
and 1.45GeV. The binding energy will increase to 99.26 MeV when the cutoff parameter is 1.60 GeV. The contribution of 
the dominant channel, A C A C ( Pq), is 90.73% ~ 81.72% when the cutoff parameter is around 1.36 ~ 1.45 GeV. The channel 
( 3 -Po) provides a fairly small contribution, less than 1%. For the corresponding bottomed state, its binding energy is 0.50 ~ 
34.63 MeV when the cutoff parameter is 0.95 ~ 1.10 GeV, which may also be a good molecule candidate. 

AgAg [/(7 P ) = 0(1)]. We consider seven channels in this case, which are listed in Table HV1 We obtain a shallow bound 
state with binding energy 0.91 ~ 20.27 MeV in the charmed sector when the cutoff parameter is 1.45 ~ 1.50 GeV. And the 
contribution of the dominant channel, A f A c ( 3 i 3 i), is 79.15% ~ 67.55%. If the cutoff parameter increases to 1.60 GeV, the 
binding energy will reach 79.78 MeV. The channel with the second largest contribution is 2 ( .E c ( 3 Pi), with a probability of 
14.55% ~ 22.90% when the cutoff parameter is around 1 .45 ~ 1.50 GeV. However, the probabilities of the other three channels, 
£ c £*( 3 Pi), E c S*( 5 Pi) and E £ .E*( 5 ,Fi), are tiny as shown in Table IXl The situation of the bottomed case is similar to that of the 
charmed case except that the binding of the former is deeper. 



B. Numerical Results for The A Q A e System with The OBE Potential 



In this subsection, we investigate the AgAg system with the OBE potential which not only includes the long-range n exchange 
interaction but also the medium-/short-range 77, cr, p and a> exchange interaction. The numerical results are shown in Table IXfl 

We obtain a weakly bound state for A C A C = 0(0 + )J. The binding energy is 2.53 ~ 55. 1 1 MeV when the cutoff parameter 

is around 0.80 ~ 1.00 GeV. Accordingly its root-mean-square radius is about 2.31 ~ 0.73 fm, which is comparable with 
the size of the deuteron. Similar to the OPE potential case, the A c A c ( 1 5o) component is dominant with a probability about 
98.69% ~ 86.79%, and the total contributions of the other channels are less than 15%. For the state A^A/, [0(0 + )], the binding 
energy is much larger as expected. Its binding energy is 27.30 MeV when the cutoff parameter is 0.80 MeV. When the cutoff 
parameter is 1.00 GeV, the binding energy reaches as high as 148.17 MeV. 

For the state A c A e = 0(0")], the binding energy is 4.69 MeV when the cutoff parameter is 1.15 GeV. When we increase 

the cutoff parameter to 1.25 GeV, the binding energy is 61.36 MeV. The probability of the dominant channel A f .A c ( 3 f ) o) is about 
85.12% ~ 71.08%. The contribution of the second dominant channel is I. C Z* C ( 3 P ), with a probability about 1 1 .48% ~ 23.57%. 
The results of the bottomed state A^Af, [0(0~)] are similar to those of the charmed case, but with deeper binding energy. The 
binding energy is 2.80 ~ 100.54 MeV with the cutoff parameter around 0.85 ~ 1.05 GeV. The probabilities of the channels 
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TABLE X: The binding solutions of AgAg with OPE potential. "A" is the cutoff parameter. "B.E." and "r rmj " are the binding energy and the 
root-mean-square radius, respectively. "P" is the probability of the individual channel which are given in Table II Vl 



KJ P ) 


Q=c 


Q=b 




A(GeV) 


0.90 


0.95 


1.00 


1.10 


1.20 


0.80 


0.90 


0.95 


1.00 


1.10 




B.E.(MeV) 


1.11 


3.99 


9.02 


26.83 


56.97 


9.09 


26.99 


40.72 


58.11 


105.17 




r„,„(fm) 


3.20 


1.89 


1.39 


0.95 


0.74 


1.08 


0.77 


0.68 


0.62 


0.52 


0(0 + ) 


Pi(%) 


98.63 


97.25 


95.69 


92.07 


88.01 


96.27 


92.77 


90.75 


88.59 


84.07 






















P 2 (%) 


0.08 


0.19 


0.35 


0.81 


1.48 


0.16 


0.42 


0.61 


0.86 


1.46 




P 3 (%) 


0.08 


0.20 


0.36 


0.87 


1.64 


0.25 


0.67 


0.98 


1.37 


2.37 




P 4 (%) 


0.86 


1.68 


2.57 


4.44 


6.29 


2.27 


4.20 


5.24 


6.28 


8.27 




P s (%) 


0.35 


0.68 


1.04 


1.81 


2.58 


1.05 


1.94 


2.42 


2.91 


3.84 




A(ueV ) 


1.36 


1.40 


1.45 


1.50 


1.60 


0.95 


1.00 


1.10 


1.20 


1.30 




B.E.(MeV) 


0.19 


9.14 


24.97 


45.32 


99.26 


0.50 


7.66 


34.63 


79.08 


142.32 






2.39 


1.06 


0.84 


0.73 


0.60 


1.94 


1.03 


0.72 


0.59 


0.50 




P x {%) 


90.73 


85.65 


81.72 


78.45 


73.05 


95.70 


91.97 


85.76 


80.01 


74.89 


0(0") 


Po(%) 


2.86 


4.24 


5.13 


5.76 


6.56 


1.26 


2.27 


3.70 


4.75 


5.46 




Pl(%) 


0.08 


0.15 


0.24 


0.35 


0.60 


0.01 


0.02 


0.10 


0.25 


0.48 




p 4 (%) 


3.97 


6.42 


8.54 


10.42 


13.74 


1.35 


2.74 


5.55 


8.59 


11.56 






2.36 


3.54 


4.37 


5.02 


6.04 


1.69 


3.00 


4.89 


6.40 


7.62 




A(GeV) 


1.45 


1.47 


1.50 


1.55 


1.60 


1.07 


1.10 


1.15 


1.20 


1.30 




B.E.(MeV) 


0.91 


7.60 


20.27 


46.97 


79.78 


1.39 


7.26 


22.28 


43.68 


105.69 




>W(fm) 


1.62 


1.99 


0.79 


0.66 


0.59 


1.40 


0.96 


0.74 


0.63 


0.51 




Pi(%) 


79.15 


73.19 


67.55 


60.70 


55.55 


91.48 


87.06 


80.52 


74.33 


63.52 


0(1") 


P 2 (%) 


14.55 


18.82 


22.90 


27.85 


31.51 


4.37 


6.89 


10.86 


14.79 


21.86 






















P 3 (%) 


2.70 


3.54 


4.40 


5.57 


6.57 


1.14 


1.81 


2.92 


4.05 


6.25 




P 4 (%) 


3.53 


4.36 


5.04 


5.77 


6.26 


2.97 


4.18 


5.63 


6.73 


8.26 




P 5 (%) 


0.08 


0.09 


0.10 


0.11 


0.11 


0.04 


0.06 


0.08 


0.09 


0.10 




Pe(%) 


0.00 


0.00 


0.00 


0.01 


0.01 


0.00 


0.00 


0.00 


0.00 


0.00 




Pi(%) 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 



A b A b ( 3 Pi) and SftS^Pi) are about 95.93% ~ 73.92% and 1.45% ~ 18.90% respectively. 

The state A C A C [I(J P ) = 0(1"")] with binding energy around 1.35 ~ 65.05 MeV and cutoff parameter around 1.16 ~ 1.25 GeV 
may also be a loosely bound state. But the binding solutions depend sensitively on the cutoff parameter. The binding energy of 
the state A h A b [0(1 ")] is 0.24 ~ 74.65 MeV when the cutoff parameter is 0.90 ~ 1 .05 GeV. 

Besides the transition induced by the OPE force in the flavor space, we have also considered the transitions caused by the 
eta meson and rho/omega meson exchange, which greatly enhances the non-diagonal matrix element in the Hamiltonian. With 
the same cutoff parameter, we can clearly see that the binding energy in the OBE case is larger than that in the OPE case. For 
example, the binding energy for the A C A C U(J P ) = 0(0 + )J state is 1.11 MeV in the OPE case if one fixes the cutoff at 0.90 GeV. 
However, the binding energy increase to 16.61 MeV in the OBE case with the same cutoff. In other words, the medium- and 
short-range attractive force plays a significant role in the formation of the loosely bound A c A f and A b A b states. 

VI. DISCUSSIONS AND CONCLUSIONS 

We have investigated the possible deuteron-like molecules composed of two heavy flavor baryons with the form of "AqBq". 
We have also performed an extensive analysis of the A e A e (Q=b,c) system, which is the heavy analogue of the H dibaryon. 

The weakly bound states are usually very sensitive to potential details including the coupling constants and form factors etc. 
Sometimes small change of the coupling constants may dismantle the bound state. 

Throughout this work, we have adopted the root-mean-square radius and binding energy of the system to judge whether the 
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TABLE XI: The binding solutions of AqAq with the OBE potential. "A" is the cutoff parameter. "B.E." and "r rm " are the binding energy and 
the root-mean-square radius, respectively. "P" is the probability of the individual channel which are given in Table llVl 



I{J P ) 


Q=c 


Q=b 




A(GeV) 


0.80 


0.85 


0.90 


0.95 


1.00 


0.80 


0.85 


0.90 


0.95 


1.00 




B.E.(MeV) 


2.53 


8.04 


16.61 


30.18 


55.11 


27.30 


45.94 


69.61 


102.62 


148.17 




r„„,(fm) 

rUlS \ / 


2.31 


1.48 


1.13 


0.91 


0.73 


0.78 


0.66 


0.58 


0.50 


0.44 


0(0 + ) 


Pi(%) 


98.69 


97.38 


95.43 


92.26 


86.79 


95.26 


92.82 


89.32 


84.43 


78.82 






















P 2 (%) 


0.08 


0.34 


1.15 


3.01 


6.04 


0.27 


0.81 


2.24 


4.81 


8.21 




P 3 (%) 


0.06 


0.18 


0.47 


1.03 


2.68 


0.31 


0.88 


1.88 


3.51 


5.60 




P 4 (%) 


0.84 


1.51 


2.15 


2.73 


3.12 


2.89 


3.73 


4.44 


4.85 


4.84 




P 5 (%) 


0.33 


0.58 


0.81 


0.97 


1.37 


1.27 


1.58 


2.13 


2.41 


2.53 




A(ueV ) 


1.15 


1.17 


1.20 


1.23 


1.25 


0.85 


0.90 


0.95 


1.00 


1.05 




B.E.(MeV) 


4.69 


13.34 


27.65 


46.59 


61.36 


2.80 


17.97 


31.35 


59.76 


100.54 




r r ,„ * (f m) 


1.25 


0.94 


0.80 


0.70 


0.66 


1.38 


0.95 


0.76 


0.62 


0.53 




Pi(%) 


85.12 


81.38 


76.85 


73.16 


71.08 


95.93 


92.31 


87.54 


80.50 


73.92 


0(0") 


Po(%) 


0.64 


0.57 


0.37 


0.19 


0.18 


1.08 


1.56 


1.73 


1.58 


1.17 




P 3 (%) 


0.75 


1.03 


1.44 


1.84 


2.11 


0.00 


0.02 


0.11 


0.76 


1.38 




P 4 (%) 


11.48 


14.67 


18.64 


21.87 


23.57 


1.45 


3.63 


7.34 


12.99 


18.90 




Ps(%) 


2.01 


2.36 


2.70 


2.93 


3.06 


1.53 


2.48 


3.27 


4.17 


4.62 




A(GeV) 


1.16 


1.18 


1.20 


1.23 


1.25 


0.90 


0.95 


0.97 


1.00 


1.05 




B.E.(MeV) 


1.35 


11.32 


24.10 


47.32 


65.05 


0.24 


11.12 


18.82 


34.83 


74.65 




r rms {im) 


1.57 


0.94 


0.78 


0.67 


0.62 


2.28 


0.92 


0.80 


0.68 


0.55 




Pi(%) 


82.41 


75.18 


70.21 


64.57 


61.64 


96.65 


90.05 


86.85 


81.50 


72.55 


0(1") 


P 2 (%) 


9.51 


13.39 


15.96 


18.67 


19.93 


1.51 


4.71 


6.26 


8.76 


12.56 






















P 3 (%) 


1.96 


2.70 


3.17 


3.68 


3.93 


0.47 


1.56 


2.10 


2.95 


4.24 




Pa(%) 


1.77 


2.24 


2.46 


2.61 


2.65 


1.29 


2.89 


3.40 


4.03 


4.64 




P 5 (%) 


4.26 


6.36 


8.02 


10.24 


11.60 


0.08 


0.77 


1.37 


2.69 


5.85 




P 6 (%) 


0.08 


0.12 


0.15 


0.02 


0.22 


0.00 


0.01 


0.03 


0.05 


0.13 




Pi(%) 


0.01 


0.02 


0.03 


0.03 


0.04 


0.00 


0.00 


0.01 


0.01 


0.03 



system is a loosely bound molecular state. There exists another intuitive approach. The relative momentum of the loosely bound 
system p ~ -\j2jJE probes distance around -, where /i is the reduced mass and E is the binding energy. For a loosely bound 
state, i should be much larger than the interaction range of the potential, which is around (0.2 - 0.3) fm. In other words, 

the size of the system should be larger than (0.6 ~ 1.0) fm. Accordingly, the binding energy should be much smaller than 
where Mb is the charmed or bottomed baryon mass. Numerically, the binding energy should be much less than 240 MeV and 
100 MeV for the charmed and bottomed systems respectively. In other words, those states in Tables IVlIllXIl which do not satisfy 
the above criteria should not be regarded as the loosely bound molecular states. 

For the spin-singlet systems with the "Agfig" form, the pion exchange force is not strong enough to form bound states for 
the five systems, [SgHg]^^ 2 ' 3 ^ 2 ', [LqAq]' s Z} q and [EqEq]^Zq' 1 ^ (Q=b,c). When we add the contributions from the scalar and 
vector meson exchanges, some bound states appear. The following five states [I^Qc]^, [EcAJ^l'g 2 , [EbAt,] 1 ^ 2 , [ScEf.]^' 2 
and [A c Q f ]£l are all very loosely bound with small binding energies and large root-mean-square radii with cutoff parameter 
1.0 ~ 1.50 GeV. They are good candidates of molecules. 

' 7=1/2 

In the spin-triplet case, the numerical results with the one-pion-exchange potential alone indicate that [EgHg] s=1 / (Q=b,c), 
[2()A2]^ 1 1 (Q=b,c), [XfrH^]^! 3 ^ 2 and [E^E^]^ may be loosely bound states. They have shallow binding solutions when the 
cutoff parameter is around 0.90 ~ 1 .50 GeV. The three states [EcS^] 1 ^ 2 , [EcE^]^ and [E&E^J^l 1 ] do not have binding solutions 
until the cutoff parameter reaches 1.80 GeV. When taking the vector and scalar boson exchanges into account, the numerical 

' 7=3/2 

results do not change significantly except that the bound state of [S C E C ] S=1 disappear and the binding energy in some channels 
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becomes deeper. Therefore, we conclude that the long-range one-pion-exchange interaction plays an dominant role in forming 
these bound states. Comparing the results of OPE model with those of the OBE model, one notices that the contribution of the D 
wave is smaller for the latter, which implies that the S-D mixing mainly comes from the pion exchange. Our results suggest that 
the states [H' G Ag]£iy 2 (Q=b,c) and [H G Q G ]£iy 2 (Q=b,c) with shallow binding solutions and reasonable cutoff parameter may 
also be good candidates of molecules. 

For the heavy analogue of the H dibaryon, our results indicate that A@ A@(Q=b,c) with quantum numbers I{J P ) = 0(0 + ), 0((T) 
and 0(1~) may all be molecules. The binding solutions of AqAq system with the OPE potential mainly come from the coupled- 
channel effect. Besides the transition induced by the OPE force in the flavor space, we have also considered the transitions 
caused by the eta meson and rho/omega meson exchange. With the same cutoff parameter, the binding energy in the OBE case 
is larger than that in the OPE case. The medium- and short-range attractive force plays a significant role in the formation of the 
loosely bound A c A f and A^A/, states. 

The authors studied the AqAq system at the quark level, and obtained bound states with mass 4516 MeV for A C A C and 9175 
MeV for A/, A;, l30ll . Theoretical investigations of these molecular states with other phenomenological models are desirable. 

If these states really exist as molecules, once produced, they will be very stable because this system decays via weak inter- 
action. It is difficult to produce the states with double charm or double bottom experimentally. However, there is still hope to 
search for these interesting long-lived molecular states with double heavy flavor at facilities such as the Large Hadron Collider 
and RHIC. 

All the molecule states (except those with Eg) are very stable because their components have a long lifetime around 10~ 13 ~ 
10~ 12 s. On the other hand, the width of S c is about 2.2 MeV [28], this narrow width ensures relatively long lifetime for the "E<JT'- 
type molecules. Such states can decay into XA+n followed by A+ — > pK~n + . For the bottomed case, I, b X -> XA° b n followed by 
A" — > A^7r~ and A+ — > pK~n + . These decay modes may be helpful to search for such states in the future experiment. 
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VII. APPENDIX 



A. Some Helpful Functions 



The functions H, etc are defined as, 



A rB 2 A rA 2 B 2 
H (A, m, r) = Y{ur) - -Y(Ar) - — Y(Ar), H^A, m, r) = Y{ur) - -Y(Ar) ^-Y(Ar) 

^ ~ ■ — A? Aj3 2 



H 2 (A,m,r) = Z x (ur) - -r Z, {Ar) - Y{Ar), H 3 (A,m, r) = Z(ur) - —Z(Ar) - -?-rZ 2 (Ar), 
M J 2 / u- i m j 2h j 



M (A, m,r) = -1 [cos(Or) - e~ Ar ] + Mi (A, m, r) = -1 [cos(0r) - «~*] - ^e~ Ar , 



M^{A,m, r) 



Or 
cos (9r) 



£ 

28A 

3 sin (Or) 3 cos (Or) 



Or 



1 A 3 AB 

T Z(Ar) - -!-rZ 7 (Ar), 

Or 3 20 3 



where, 



and 



P 2 - A 2 - m 2 , u 2 = m 2 - Q 2 , 2 = -(m 2 - Q 2 ), A 2 = A 2 - Q 2 , 



Y(x) = —, Z(jc) = |1 + - + ^-|7(x), Z l (x) = \- + \\Y(x), Z 2 (x) = (1 + x)Y(x) 



Fourier transformation formulae read: 



1 



u 2 +Q 2 

Q 



u 2 + Q 2 4tt 



u Q 2 u 

—H (A, m, r), -> -— H,(A, m, r), 

4-n u + Q 2 4jt 

iu 3 „.. , QiQj 
— rH 2 (A, m, r), 



u 2 + Q 2 Yin 



[H 3 (A, m, r)kij + ffi(A, m, r)6ij\ , 



(26) 



(27) 



where, k,j = 3-V- - 5y. If the form factor is not introduced, there will be delta terms in the second and fourth formule. In the 

above expressions, we employ another function, -^-[H\(A, m, r) - H (A, m, r)], to substitute for the delta term. If one neglects 
the delta term, one should adopt the following formulae, 



Q- 



u + Q 2 



~—H (A,m,r), 

4-7T u 2 +Q 



TZTa ~* — jy- [Hs(A, m, r)k u + H (A, m, r)6ij\ . 



If u 2 = m 2 x - Qg < 0, the last formula of the Eq. (|27T > should be 



QiQj 



r i 

-— [M 3 (A, m, r)kij + Mi (A, m, r)5 u \ . 



u 2 + Q 2 I2n 

Accordingly, we make the replacement M\ (A, m, r) — > Mq(A, m, r) to neglect the delta term. 



(28) 



(29) 



B. The Numerical Results of The "A Q B Q " Systems When The Contact Term Is Included 



For comparison, we collect the numerical results of the "AqBq" systems in Tables IXTTl and IXUTI when the contact term is 
included. 
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TABLE XII: The binding solutions of the spin-singlet "AgBg" systems when the interaction potential includes the contact term. 



5=0 



States 



A(GeV) 



Q = c 
B.E.(MeV) 



,(fm) 



A(GeV) 



Q = b 
B.E.(MeV) 



s (fm) 



[Z e H e ]j - 



[Z e H e ] s 7 



[s e n e ]f 



[H e A e ] s 



[1qSq] s 3 



y_3 



0.80 
0.85 
0.90 

1.30 
1.40 
1.50 

1.10 
1.20 
1.30 

0.90 
0.95 
1.00 

1.00 
1.10 
1.20 

1.80 
1.90 
2.00 

0.90 
0.95 
1.00 

1.50 
1.60 
1.70 



52.87 
82.54 
137.34 

0.23 
1.48 
4.11 

0.94 
1.03 
2.53 

2.55 
8.80 
18.69 

0.73 
1.77 
2.31 

0.21 
0.83 
2.00 

2.59 
10.29 
24.38 

0.42 
1.28 
2.66 



0.64 
0.55 
0.45 

5.44 
3.03 
2.10 

4.56 
3.29 
2.35 

2.08 
1.26 
0.94 

3.67 
2.55 
2.28 

5.96 
3.99 
2.69 

2.12 
1.22 
0.88 

4.55 
3.06 
2.30 



0.80 
0.85 
0.90 

1.00 
1.10 
1.20 

0.90 
1.00 
1.10 

0.80 
0.85 
0.90 

0.80 
0.85 
0.90 

0.80 
1.00 
1.20 

0.80 
0.85 
0.90 

0.85 
1.00 
1.20 



162.52 
219.83 
321.82 

0.45 
0.45 
4.61 

1.12 
5.36 
7.59 

4.89 
15.30 
31.25 

0.55 
4.93 
10.39 

0.13 
0.19 
0.32 

6.13 
14.88 
30.65 

0.85 
2.25 
3.69 



0.34 
0.30 
0.26 

3.49 
3.49 
1.66 

2.39 
1.41 
1.32 

1.18 
0.78 
0.61 

2.94 
1.25 
0.96 

6.53 
6.06 
5.35 

1.08 
0.81 
0.63 

2.56 
1.82 
1.61 



[Z Q Z Q ]' S 



[E Q E Q ]' S 1 



[E Q Q Q ] S 



1.10 
1.20 
1.30 

1.10 
1.20 
1.30 

1.25 
1.30 
1.40 

1.00 
1.10 
1.20 

0.90 
0.95 
1.00 



0.45 
1.02 
1.86 

1.85 
6.37 
12.28 

0.23 
0.47 
1.25 

0.52 
1.99 
3.41 

0.20 
2.11 
5.77 



4.39 
3.32 
2.70 

2.59 
1.60 
1.27 

5.31 
4.41 
3.08 

4.15 
2.50 
2.06 

5.29 
2.33 
1.55 



0.85 
1.00 
1.20 

1.00 
1.10 
1.20 

0.95 
1.00 
1.20 

0.90 
1.00 
1.20 

0.80 
0.85 
0.90 



I. 55 
7.16 

II. 56 

2.51 
11.11 

21.83 

0.80 
1.79 
5.95 

1.34 
7.82 
15.97 

1.64 
6.53 
14.20 



2.00 
1.16 
1.08 

1.65 
1.05 
0.87 

2.65 
1.97 
1.38 

2.18 
1.21 
1.06 

1.87 
1.12 
0.86 
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TABLE XIII: The binding solutions of spin-triplet "AqBq" systems when the interaction potential includes the contact term. 



States 


A(GeV) 


B.E.(MeV) 


Q = c 


P S (%) 


S 

Pd(%) 


= 1 

A(GeV) 


B.E.(MeV) 


Q = b 


Ps(%) 


Pd(%) 




0.90 
1.00 
1.10 


1.39 
8.64 
20.31 


2.97 
1.51 
1.14 


96.08 
93.78 
92.74 


3.92 
6.22 
7.26 


0.80 
0.90 
0.95 


5.24 
18.07 
27.54 


1.51 
1.02 
0.90 


88.01 
86.04 
85.53 


11.99 
13.96 
14.47 


, r _3 


0.85 
0.90 
0.95 


3.04 
8.66 
16.21 


1.98 
1.30 
1.00 


97.19 
96.77 
97.31 


2.81 
3.23 
2.69 


0.80 
0.85 
0.90 


16.92 
36.05 
56.13 


0.82 
0.65 
0.54 


92.42 
91.49 
92.82 


7.58 
8.51 
7.18 




0.90 
n 

1.00 


2.37 
1 o ^ 

20.87 


2.12 

1 1 Q 

1 . ly 

0.91 


99.38 

QS S 1 

yo.J 1 

98.08 


0.62 

1 zLQ 
1 .t+y 

1.92 


0.85 
n on 

0.95 


11.23 

11 QA 

J 1 .70 

56.35 


0.80 
n ai 

0.52 


99.56 

y 1 JO 

95.42 


0.44 

7 A 7 

4.58 


[SgQ^]^ -1 


0.95 
1.00 
1.10 


1.91 
5.15 
15.31 


2.42 
1.62 
1.08 


99.84 
99.71 
99.40 


0.16 
0.29 
0.60 


0.80 
0.85 
0.90 


1.35 
5.53 
12.11 


2.03 
1.20 
0.92 


99.27 
99.48 
99.17 


0.23 
0.52 
0.83 




1.00 
1 1 n 

1.20 


0.73 

1 77 

2.31 


3.67 
7 

2.28 


100.00 
i on nn 

1UU.UU 

100.00 


0.00 
u.uu 

0.00 


0.80 
n 

0.90 


0.55 

A. CH 

10.39 


2.94 

1 .ZJ 

0.96 


100.00 
1 nn nn 

100.00 


0.00 
n nn 

u.uu 
0.00 




0.80 
0.85 
0.90 


2.19 
5.32 
9.34 


2.50 
1.75 
1.39 


95.11 
95.23 
95.67 


4.98 
4.77 
4.33 


0.80 
0.85 
0.90 


28.86 
41.93 
56.53 


0.79 
0.67 
0.58 


87.95 
89.83 
92.07 


12.05 
10.17 
7.93 


_ ?= i 
[£ e H e ] s 2 


0.90 
1.00 
1.10 


0.57 
9.52 
25.78 


3.99 
1.31 
0.92 


99.96 
99.88 
99.88 


0.04 
0.12 
0.12 


0.80 
0.85 
0.90 


1.10 
6.44 
15.96 


2.22 
1.15 
0.84 


99.12 
99.72 
99.56 


0.08 
0.28 
0.44 


'=4 

[£gH e ] s 2 


0.95 
1.00 
1.10 


0.49 
1.56 
6.02 


4.16 
2.60 
1.45 


99.81 
99.68 
99.61 


0.19 
0.32 
0.39 


0.80 
0.85 
0.90 


4.04 
0.90 
14.91 


1.29 
0.98 
0.82 


99.59 
98.73 
97.99 


0.41 
1.27 
2.01 


[S e A e ]f ! 


0. 95 

1 nn 

1. UU 

1.10 


1.25 

1 A/1 

12.23 


2.91 

L.oJ 

1.14 


99.88 

yy.od 

99.83 


0.12 
n i 7 

0.17 


0.80 
n 

0.90 


3.15 

C /in 
o.4U 

15.51 


1.43 

1 ni 
1.U1 

0.82 


99.76 
qq /in 

99.15 


0.24 
n An 

0.85 


[H e H Q ]J =0 


0.90 
n 

1.00 


4.89 
1 a An 

28.43 


1.67 
i i n 

0.86 


98.42 

70.7J 

99.36 


1.58 

1 07 

1 .U / 

0.64 


0.80 

U.oJ 

0.90 


5.71 

1 A A7 

35.81 


1.39 
n on 

0.66 


86.06 

Q1 1A 

96.05 


13.94 

S 7A 

o.zo 
3.95 


[S e H e ]J =1 


0.90 
0.95 
1.00 


0.81 
4.39 
9.88 


3.41 
1.69 
1.22 


99.63 
99.39 
99.31 


0.37 
0.61 
0.69 


0.80 
0.85 
0.90 


0.77 
10.07 

23.63 


2.50 
0.93 
0.70 


98.86 
98.57 
98.19 


1.14 
1.43 
1.81 




0.90 
n 

1.00 


0.45 

^ QA 

15.35 


4.22 

1 AQ 

1.04 


99.75 

QQ 1A 

99.23 


0.25 
n aa 

0.77 


0.85 
n on 

0.95 


3.01 

1 Q AA 
17.44 

41.92 


1.45 

n 7A 
u. /o 

0.58 


97.79 
Q7 nn 

y 1 .UU 

97.13 


2.21 
i nn 

2.87 


[A fl Q fl ]5 =0 


0.90 
0.95 
1.00 


0.20 
2.11 
5.77 


5.29 
2.33 
1.55 


100.00 
100.00 
100.00 


0.00 
0.00 
0.00 


0.80 
0.85 
0.90 


1.64 
6.53 
14.20 


1.87 
1.12 
0.86 


100.00 
100.00 
100.00 


0.00 
0.00 
0.00 
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C. The Numerical Results of The Baryon-antibaryon Systems With the OPE potential 



As a byproduct, we present the binding solutions of the heavy baryon-antibaryon systems with the pion-exchange potential. It 
is straightforward to obtain the potential via changing the sign of the potential for the baryon-baryon systems since the G-parity 
of the pion is negative. Our results indicate that the pion-exchange alone is strong enough to form some bound states. The 
numerical results are collected in Tables IXIVUXVl 



TABLE XIV: The numerical results of the spin-singlet heavy baryon-antibaryon systems with the OPE potential. "A" is the cutoff parameter. 
"B.E." is the binding energy, and "r rmi "is the root-mean-square radius which reflects the size of the bound state, "x" denotes no binding 
solutions 



Systems 



A(GeV) 



Q=c 
B.E.(MeV) 



S=0 



A(GeV) 



Q=b 
B.E.(MeV) 



r rms (fm) 



1.80 
2.00 
2.50 



0.01 
0.02 
0.03 



6.75 
6.66 
6.50 



ry 5' i'= 3 / 2 



0.80 
1.00 
1.30 



0.06 
0.13 
0.29 



9.34 
4.94 
2.44 



T7T & l'=0 



D. The Dependence of The Binding Energy on The Cutoff Parameter 
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TABLE XV: The numerical results of the spin-triplet heavy baryon-antibaryon systems with the OPE potential. "A" is the cutoff parameter. 
"B.E." is the binding energy while "r„„ s " is the root-mean-square radius. "Pg" and "P " are the probabilities of the S wave and the D wave, 
respectively. "X" denotes no binding solutions. 



S=l 

Q=c Q=b 



Systems 


A(GeV) 


B.E.(MeV) 


r„,„(fm) 


Ps 


Pd 


A(GeV) 


B.E.(MeV) 




Ps(%) 


Pd(%) 


ry p< n /=l/2 

L^e^gls 


1.20 
1.30 
1.40 


I. 37 
4.95 

II. 24 


2.80 
1.69 
1.25 


89.79 
84.74 
81.80 


10.21 
15.26 
18.92 


0.80 
0.90 
1.00 


4.05 
9.86 
19.09 


1.60 
1.21 
0.98 


74.87 
70.37 
67.15 


25.13 
29.63 
32.85 



ry ,1=3/2 






X 






1.20 
1.30 
1.50 


0.72 
1.75 
5.80 


2.68 
1.88 
1.18 


94.08 
92.19 
89.44 


5.92 
7.81 
10.56 




1.60 
1.80 
2.00 


0.03 
1.38 
5.87 


6.65 
2.99 
1.56 


97.47 
94.34 
92.02 


2.53 
5.66 
7.98 


0.90 
1.10 
1.30 


I. 41 
4.42 

II. 71 


2.75 
1.59 
1.04 


87.92 
86.20 
84.88 


12.08 
13.80 
15.12 




2.10 
2.30 
2.50 


0.46 
3.16 
9.15 


4.13 
2.89 
1.20 


97.15 
94.51 
92.55 


2.85 
5.49 
7.45 


1.00 
1.20 
1.40 


0.24 
1.89 
6.22 


4.52 
1.91 
1.20 


93.44 
90.07 
87.85 


6.56 
9.93 
12.15 


[E Q E Q ]'f' 






X 










X 








FIG. 5: The dependence of the binding energy on the cutoff parameter for the "AqBq" system with the OPE potential. 
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A I GeV) A (GeV) A (GeV) A (GeV) 

FIG. 6: The dependence of the binding energy on the cutoff momentum for the spin-singlet system "A C B" with the OBE potential. 




_i i i i_ 



0.92 0.94 0.96 0.98 1 
Symbol L (GeV) 



0.8 0.85 0.9 0.95 1 1.05 1.1 
A (GeV) 



FIG. 7: The dependence of the binding energy on the cutoff parameter for the spin-singlet "A^Bfy" system with OBE potential. 
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A(GeV) A (GeV) A (GeV) A(GeV) 

FIG. 8: The dependence of the binding energy on the cutoff parameter for the spin-triplet "A C B C " system with the OBE potential. 




Q F ■ ... I .... I ■ ... I ■ ... I .... 1 \ I I I I I I I 1 E. ... I ■ ... I .... I ■ ... I .... I . ... : 

0.8 0.B5 0.9 0.95 1 1.05 0.9 0.95 1 1.05 1.1 1.15 1.2 1.25 1.3 0.9 1 1.1 1.2 1.3 1.4 1.5 

A (GeV) A (GeV) A (GeV] 



FIG. 9: The dependence of the binding energy on the cutoff parameter for the spin-triplet u A b B b " system with the OBE potential. 
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FIG. 10: The dependence of the binding energy for A^A^ on the cutoff parameter with the OPE potential. 



